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Abstract 

In this paper we examine the wave properties of hot plasma liv- 
ing in Schwarzschild magnetosphere. The 3+1 GRMHD perturba- 
tion equations are formulated for this scenario. These equations are 
Fourier analyzed and then solved numerically to obtain the disper- 
sion relations for non-rotating, rotating non-magnetized and rotating 
magnetized plasma. The wave vector is evaluated which is used to cal- 
culate refractive index. These quantities are shown in graphs which 
are helpful to discuss the dispersive properties of the medium near the 
event horizon. 
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1 Introduction 

In general relativity (GR), the term black hole evokes the mysterious gravity. 
Black holes are among the most remarkable predictions of GR, which are a 
reality today. These are usually found in X-ray binaries and in the centers of 
galaxies. The existence of rotating stars indicates their rotational behavior. 
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A rotating black hole gradually reduces to a Schwarzschild (non-rotating) 
black hole by the extraction of its rotating energy. Plasmas are abundant 
in nature. More than 99% of all known matter is in the plasma state. All 
the stars are made of plasma, and even the space between the stars is filled 
with plasma. Commonly space plasma occurs in a hot state. The strong 
gravity of the black hole strips the plasma from the surrounding star. Thus 
the plasma gathered around the black hole in the form of accretion disk. The 
moving plasma creates a magnetic field. The region surrounding the black 
hole admitting the magnetic field is known as magnetosphere. The theory of 
general relativistic magnetohydrodynamics (GRMHD) is probably the most 
accurate approach to investigate the dynamics of relativistic, magnetized 
plasma. 

The Schwarzschild black hole is non-rotating and hence the magneto- 
spheric plasma falls freely along radial direction only. Perturbations in the 
Schwarzschild regime, either geometrical or physical, have always been of 
interest by the relativists. Regge and Wheeler [T], Zerilli [2] and Price [3] 
discussed the gravitational perturbations. Fiziev [I] presented the exact solu- 
tion of Regge- Wheeler equations. These equations describe the axial pertur- 
bations of the Schwarzschild metric in linear approximation. The quasi-static 
problem of electric field in the Schwarzschild black hole was solved by Hanni 
and Ruffini [5J. Sakai and Kawata [6] developed a special relativistic ap- 
proach for a linearized treatment of plasma waves in the Schwarzschild black 
hole magnetosphere. 

In general relativity, a 3+1 hypersurface split of spacetime is appropriate 
to understand the black hole physics. This split was developed by Arnowitt, 
Deser and Misner (ADM) [7] to study the quantization of gravitational field. 
The formalism has wide applications in numerical relativity. The 3+1 ap- 
proach has been used by many authors [S]-[in] to discuss different features in 
GR. Thorne and Macdonald [H]-[l2] developed the electromagnetic theory in 
black hole regime using this formalism. Holcomb and Tajima [13], Holcomb 
[H] and Dettmann et al. [15] studied some properties of wave propagation 
for the Friedmann universe. Buzzi et al. [H] investigated relativistic two 
fluid plasma wave properties in the vicinity of the Schwarzschild black hole. 
Ali and Rahman [T7] adopted the technique used by Buzzi et al. [16] to 
analyze the transverse electromagnetic waves propagating in a plasma close 
to the Schwarzschild-de Sitter black hole. Zhang [18] formulated the black 
hole theory for stationary symmetric GRMHD. He [19] also discussed the 
behavior of cold plasma perturbations in the Kerr magnetosphere. Sharif et 
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al. [2U]-[24J discussed properties of plasma waves by using real and complex 
wave numbers. The analysis was given both for cold and isothermal plasmas. 

A lot of work has been done using cold and isothermal plasma but no one 
has used hot plasma in this context which is the basic constituent of nature. 
This is the most general plasma which reduces to cold and isothermal plasma 
with some restrictions. We have considered this plasma around the black hole 
to check the possibility of receiving information. This work focusses on the 
investigation of hot plasma wave properties in the Schwarzschild magneto- 
sphere. We shall apply perturbation and Fourier analysis techniques. The 
dispersion relations are calculated with the help of the software Mathematica 
to obtain the wave vector. This will be used to evaluate the refractive index 
and its change with respect to angular frequency. The wave properties will 
be found through these quantities. 

The outline of the paper is as follows. In Section 2, we shall provide the 
general line element and restrict it to the Schwarzschild planar analogue. Sec- 
tion 3 includes the plasma assumptions for perturbation and Fourier analysis. 
Moreover, the perturbed and Fourier analyzed 3+1 GRMHD equations are 
specified for hot plasma. In Sections 4, 5 and 6, we restrict these equations 
to non-rotating, rotating non-magnetized and rotating magnetized plasmas 
and discuss the wave properties. Section 7 contains summary of the results. 



2 3+1 Split of Spacetime 

In ADM 3+1 split, the four dimensional spacetime is decomposed into a 
succession of three dimensional spacelike hypersurfaces with the directions 
normal to them taken to be universal time direction. The line element of 
spacetime in 3+1 formalism can be written as [Tj5] 

ds 2 = -a 2 dt 2 + 7y(c^ + (3 i dt)(dx j + (3 j dt), (2.1) 

where a is the lapse function which describes the ratio of fiducial proper 
time to universal time, i.e., a = 4^. The shift vector components j3 % indicate 
the shift of spatial coordinates as one moves from one hypersurface to next. 
jij are components of three dimensional hypersurface (the absolute space) 
metric. A natural observer is associated with this spacetime, called the 
fiducial observer (FIDO). 

The Schwarzschild black hole is non-rotating, thus in the planar analogue 
of Schwarzschild geometry, the shift f3 vanishes and the above line element 
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becomes 

d s 2 = -a 2 (z)dt 2 + dx 2 + dy 2 + dz 2 . (2.2) 

Here, the directions z, x and y are analogous to Schwarzschild's r, <fr and 9 
respectively. 

3 3+1 GRMHD Equations with Relative As- 
sumptions 

For the plasma existing in the general and Scwarzschild planar analogues 
(given by Eqs. (l2.1l) and (12.21) ). the 3+1 GRMHD equations are given in 
Appendix A. We assume that hot plasma surrounds the Schwarzschild black 
hole. The specific enthalpy of the fluid is 



P + P fOl\ 

fi = , 3.1 

Po 

where po, p and p denote the rest-mass density, moving mass density and 
pressure respectively. Equation (13. ip indicates the exchange of heat between 
the plasma and the magnetic field of the fluid. 

We can modify the 3+1 GRMHD equations given by Eqs. (lA6j) - (lA10l) for 
hot plasma living in the Schwarzschild spacetime as follows 

_ = Vx(aVxB), (3.2) 
V.B = 0, (3.3) 

%£±$- + (p + P)7 2 V.^ + ( P + p^V.iaV. V) V 

+ (p + p)V.(aV) = 0, (3.4) 

> + pW + %) k + (p + ph% Vj - ± Bi B 3 } g* 



+V.V) V + 7 V,(V.V)(p +p)-[ ^8 i3 - -^—BiBj ) V{V k 



4tt j 4tt 

-(p + ph 2 ai - Pi . + J-(V x B)iV.(V x B) - -1j(qB) 
Ait Qixa z 



+^(aB i ) J Bi - -L[B x {V x (V x (aV x B))}],, (3.5) 
\^{{P + Ph 2 ~P} + (P + p)7 2 V.a + (p + p) 7 4 V.(V.V) V 



_i_ a .{(V x B) x B} = i-B.{V x (V x B)} 

<^ 2-(V x B). — (V x B) + (V x B).(a x B) 

An y a dt ' 

where B and V are the velocity and magnetic field of the fluid as measured 
by the FIDO. 

We consider rotating background (plasma is rotating), plasma is not only 
moving along radial direction but it rotates as well. Due to its rotation, wave 
propagates along axial direction also and thus propagates in (x, z)-plane. It 
is assumed that the FIDO measured velocity of fluid and magnetic field lie 
in xz-pl&ne 

V = V(z)e x + u(z)e z , B = B[X(z)e x + e z ), (3.7) 
where B is an arbitrary constant. Here A, u and V are related by [20] 



V F 

V = + Xu, (3.8) 

a 

where V F is an integration constant. The Lorentz factor 7 = , 1 takes 
the following form 

The plasma flow in the magnetosphere can be characterized by its density 
p, pressure p, velocity V and magnetic field B. When the flow is perturbed, 
these variables will become 

p = p° + 5p = p° + pp, p = p° + Sp = p° + pp, 

V = V° + 5V = V° + v, B = B° + 5B = B° + J Bb, (3.10) 

where p°, p°, V° and B° are unperturbed quantities. The linear perturba- 
tions in these quantities are denoted by 5p, 8p, 5~V and <5B respectively. 

We shall use the following dimensionless notations for the perturbed quan- 
tities 

p = p(t,z), p = p(t,z), 

v = SV = v x (t, z)e x + v z (t, z)e z , 

b = — = b x (t,z)e x + b z (t,z)e z , (3.11) 



where p, p, v x ,v z , b x and b z are dimensionless quantities. For Fourier analysis, 
we assume the harmonic space and time dependence of perturbations 



p(t, z) 


= C\t 


-L(u)t- 


-kz) 


p(t, z) 


= c 2 e-^ l - kz \ 




v z (t,z) 


= c 3 e~ 




-fcz) 


v x (t,z) 


i = c 4 e- L{ut - kz) , 




b z (t,z) 


= c 5 e" 


-i(yit- 


-fcz) 


b x (t,z) 




(3.12) 



where k is z-component of the wave vector (0, 0, k) and u is the angular 
frequency. We can use the wave vector to obtain the refractive index which 
we shall use to deduce the plasma wave properties. We define the wave vector 
and the refractive index as follows: 

• Wave Vector: A wave vector is a vector which points in the direction 
of propagation of wave. Its magnitude gives the wave number. 

• Refractive Index: It is the ratio of the speed of light in vacuum to 
the speed of light in the other medium [25]. If the refractive index is 
greater than one and its change with respect to angular frequency is 
positive, the dispersion is said to be normal, otherwise anomalous. 

When we introduce linear perturbations from Eq. fl3.10p . the perfect GRMHD 
Eqs.(l3T2l- (j3TB|l turn out to be 

?2?2 = V x (cvv x B) + V x (qV x SB), (3.13) 
at 

V.(5B)=0, (3.14) 

I^p + M + + 19 + + + 

a at a at 

= -2(p + p) 7 2 (V.v)(V.V)ln 7 - (p + p) 7 2 (V.VV).v 

+(p + p)(v.Vlnw), (3.15) 

B 2 \ 111 dv j 

(p + P)7 2 + ^ J 5, + (p + pHV^ - -B lBj } - - 

+i-[B x {V x + (P + ptfvijV* + (p + p) 

x-f 4 ViV jjk V j V k + 7 V l (V.V)(5p + 8p) + 7 2 ^(v.V)(p + p) 

+ 7 2 ^(V.V)(p + p) + 7 4 (2V.v)^(V.V)(p + p) - -L^aSB^ 

Ana 
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-7 2 {(5p + 5p)V ] + 2(p + p) 7 2 (V.v)^' + (p + pyyVij 
- 7 4 ^{(<Sp + 5p)V j + 4(p + p) 7 2 (V.v)^' 

+ (p + P y}V 3:k V k , (3.16) 

f i dispel + 2_ {p+ * -im +{6p+ 5pW 

a at a at a at 

x [(V.a) + 7 2 V.(V.V)V] + (p + p) 7 2 [v.a + 2 7 2 (V.v)(V.a) 

+4 7 4 (V.v)V.(V.V)V + 7 2 (v.(V.V)V + V.(v.V)V + V.(V.V)v)] 
= — [cSB.V x (V x B) + B.V x (v x B) + B.V x (V x SB) 

(v x B + V x <5B).— (V x B) - -(V x B). ( x B 

+V x ^5 j _ ( v x B + V x 5B). (a x B) - (V x B).(a x (SB) 



9* / 

5B dB B d(5B) 1 



a.[(v x B) x B + (V x 5B) x B 



a at a' at 2ir 
+ (V x B) x <SB]. (3.17) 

When we use assumptions given by Eq. fl3.lll) . the component form of 
Eqs. fl3.13p - fj3.17n can be written as 

1 db x 

— 7tt + ub x z = (v x - Xv z + Vb z - ub x )V In a 
a at 

+K, - Xv z , z - X'v z + V'b z + Vb z>z - u'b x ), (3.18) 
X db 

-^ + ub ZtZ = 0, (3.19) 
a at 

b z , z = 0, (3.20) 

1 dp 1 dp . . ~ r .\dv x . , . 2 

P-T^7 + P-T^r + [P + P)l V{ — — + uv xz ) + (p + p)7 u 
a at a at a at 

1 ov 

X a~8t +(P + P^ 1 + 7 2 « 2 K* = ~iMp + 2 7 V 2 )K' 

2_.2\/ 1 i .,2.A U 



+2 7 WuV + (p + p)[(l - 2 7 V)(1 + 7 V) 



27VWK, (3.21) 
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(p + p)y(i + y^) + _|--^ + |(p +P )Av 4ff 

x ^ + { (p + p) 7 2 (l + 7V 2 ) - ^ } + { (p + p) 7 W 

+— ; — 1 uv z , ; — (1 — u 2 )b x , — \a'(l — u 2 ) — auu'\ b x 

An J z ' z 4tt v ; ' Ana 1 v ; J 

+ 1 2 u(pp + pp) {(1 + 7V 2 ) V + 7 2 wW} + -f 2 uV(p'p + pp' 

+p'p + PP') + [(P + P)7 4 « { (1 + 47V 2 ) uu' + AW {I + 7V 2 ) } 

+ + 7 Vl + 2 7 V 2 ) (p' + p')R + [(p + p) 7 2 { (1 + 2-yV) 

(1 + 2 1 2 V 2 )V> - 7 W + 2 7 2 (1 + 2 7 V)«W} + ^(Aa)' 
+ 7 V(1 + 2 7 V)(p' + p')]v g = 0, (3.22) 

{ {p+p)lHl +lV) + ^} + {(,+„) - 

x ^ ^ + j (p + p)7 2 (l + 7 2 m 2 ) - j ut; Zl z + { (p + p)7 4 ^^ 

A.B 2 1 XB 2 2 Z? 2 r . . . 

H — : — f ™ I]Z + — — (1 - u )b x>z + \{oiX) + a X- uX(ua 

An J An Ana 

+u'a)} b x + (pp + pp)7 2 {a z + uu'{\ + 7V) + 7V W} 

+ (1 + 7 V)(p'p + pp') + 7 V(P'P + P~p') + [(p + P)7 4 

x{«V(l + 4 7 V 2 ) + 2F(a 2 + W(l + 2 7 V))} + 

Ana 

+2 7 V\/(p' + p')R + [(P + P)7 2 {«'(! + 7V) (1 + 47V) 

+2« 7 2 (a 2 + (1 + 2 7 VW')I - ^-%A)' + 2 7 2 «( 1 

Ana 

+ 7 V)(p' + p> 2 = 0, (3.23) 

2 [. , 4 A5 2 , "1 «% 2 5 2 rx 

+2 M (A« - V)} ^ + — {1 - 2V(Xu ~V)}^ + (pp + pp) 



x 7 2 m {a 2 + 7 2 ( W + uu') } + j (p + p) 7 4 «V 



XB 2 
An 



, . 4 2 A 2 5 2 | XB 2 u XB 2 V 
+ <{p + Pn u + — — > v Z)Z + — — b x>z b 



An J ' 1- 1- ; 

+[(P + vW + 2yW)(l + 2 7 V) + 7 W(1 + 4 7 V)} 

+^-(3Xa z + X')]v z + [(p + p)^u {V'(l + 4 7 V 2 ) + 2V 

3B 2 B 2 

(a z + 2yW)} a z X]v x + — [3a z {2Xu - V) + 2u'X 

An ' An 

XB 2 

+uX' - V']b x - -—{3Va z + V')b z = 0. (3.24) 
An 

Using time dependence harmonic space of the perturbed variables from Eq. fl3.12p . 
we get the Fourier analyzed form of Eqs. fl3.18p - fl3.24p 

c 4 (c/ + ika) — c 3 {(aA)' + LkaX} + c 5 (aV)' — c 6 {(a«)' — luj 
+ikua} = 0, (3.25) 

c 5 ( — + ika) = 0, (3.26) 
a 

c 5 ik = 0, (3.27) 

ci(— p) + c 2 (—p) + c 3 (p + p)[— 7 2 u + (1 + -i 2 u 2 )ik 
a a a 

-(1 - 2 7 V)(1 + 7 V)- + 2 7 4 «W] + c 4 (p + p) 7 2 [( — 

u a 

+iku)V + u{l + 2 1 2 V 2 )V + 2 7 2 uW] = 0, (3.28) 
Cl [p 7 2 M{(l + 7 2 V 2 )V" + j 2 Vuu'} + 7 2 V«(p' + ukp)} 
+c 2 [p 7 2 w{(l + 7 V 2 )U' + ^Vuu'} + 7 Vu(p' + iifep)] 

+c 3 [(p + P)7 2 {(1 + 2 7 V)(1 + 2 7 V 2 )V" + (— + Lku)-i 2 Vu 

a 

_ 7 W + 2 7 2 (1 + 2j 2 u 2 )uVu'} + 7 V(1 + 2 7 V)(p' + p') 

+p^(Xa)' + ^!(- + tfcu)] + c 4 [(p + p) 7 4 u{(l + 4 7 V 2 ) 
Ana An a 

xuu + AW (I + 7 V 2 )} + (p + p) 7 2 (l + 7 V 2 )(— + iku) 

a 

o / o ox / / 1 1 B 2 ua' B 2 .luj 

+ 7 y i + 2 7 2 v 2 )( P ' + p') + — - + Lku)] 

47ra 47r a 

B 2 

+c 6 [ouu' - a'(l - m 2 ) - (1 - M 2 )iA;a] = 0, (3.29) 

Cl [p 7 2 {a 2 + (1 + 7 V)W + 7 2 wW} + 7 V(p' + ukp)\ 
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+c 2 [p 7 2 R + (1 + 7 2 «>m' + t 2 «VF'} + (1 + 7V) 

x (p' + tfcp)] + c 3 [(p + p)7 2 {(l + 7 2 w 2 ) ( h Jew) 

a 

W(l + 7 V) (1 + 4 7 V) + 2w 7 2 (a z + (1 + 2 7 V) W)} 

+2j 2 u(l + 7 V)(p' + p') - ^(Aa)' - + ifcu)] 

47ra; 47T a 

+c 4 [(p + p)7 4 {(— + + u 2 V'(l + 4 7 V 2 ) + 2V(a z 

a 

+ (1 + 2 7 V)W)} + 2 7 VVV + P') + ^-{— + t Jfe u ) 

XB 2 ua' B 2 

H 1 + cJ {(AaV + a'A - u\(ua' + u'a)} 

Ana L Ana 

+^-(l-u 2 )ik] = 0, (3.30) 

47T 

Cl p7 2 [-— + u{a z + 7 2 (\/\/' + W)}] + c 3 [- — {(p + p)7 4 « 
a a 

_B 2 A A 2 S 2 
+ -— (uA - V)} + Lk{(p + p) 7 V + — — } + (p + p)7 2 {(a 2 
An An 

+2 7 W)(1 + 2 7 V) + 7 W(1 + 4 7 V)} + (3Aa z + A')^] 

47T 

Ik ) 

+c 4 [ {(p + p) 7 4 ^ - — (uA - V)} + ^{(p + p) 7 4 ^ 

a 47r 

AS 2 

' ' } + (p + p)7 4 «{^'(l + 4 7 V 2 ) + 27(a 2 + 2 7 W)} 



47T 

35 2 £ 2 r 



a 2 A] + c 6 — [ {A + 2w(uA - 7)} + tkuX + 2u'\ 



An An a 

+3a z (2Xu -V) + tiA' - V] + c 2 p[- — (7 2 - 1) + u-f{a z 

a 

W{W + uu')}} =0. (3.31) 



Equations (13.261) and (13.271) give c 5 = indicating that there are no pertur- 
bations in the ^-component of magnetic field. 

4 Non-Rotating Plasma Flow 

In non-rotating plasma flow, the magnetospheric perturbations are only along 
z-axis. The FIDO measured magnetic field and velocity admit only z- 
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component i.e., V = = A which gives v x = = b x . The Fourier analyzed 
perturbed GRMHD equations for non-rotating plasma can be obtained by 
substituting these assumptions along with the vanishing respective Fourier 
constants C4 and cq in Eqs. fl3.25l) - fl3.31l) 

--c 5 = 0, (4.1) 

a 

ikc 5 = 0, (4.2) 
ci{ p) + c 2 { p\ + c 3 {p + p){(l + 7 u 7 m 

-(1-2 7 V)(1 + 7 V)-} =0, (4.3) 

u 

ci{fri 2 {a z + W(l + 7 2 m 2 )) + 7 2 m 2 (p' + ikp)} 
+c 2 {jri 2 {a z + W(l + 7V)) + (1 + 7 2 « 2 )(P' + t*p)} 

+C 3 7 2 [(P + P){(1 + 7 2 « 2 )(— — + th)+«'(l+7V) 

x(l + 47V) + 27 2 Ma,} + 2M(l + 7 V)(p , +p / )] =0, (4.4) 

Or bid . 2 / \ ~i r / 2 -i \ 

c iP7 I V u i a z + 7 ««)} + C 2P{ (7 - 1) 

a a 

+7 2 w(a 2 + 7 2 W)} + c 3 (p + ph 2 { — 7 2 w 

a 

+7 2 u 2 ik + (a 2 + 2 7 W) (1 + 2 7 V) } = 0. (4.5) 
We shall use these equations to obtain dispersion relations. 

4.1 Numerical Solutions 

To obtain numerical solutions, we use the following assumptions 
• Time lapse: a = tanh(10z)/10, 



Specific enthalpy: p = y 1 ( tanh (^ 02 0/ 10 ) ; 

Stationary fluid: 07 = 1 =>■ 7 = 1/ y/1 — u 2 = 1/a 

For the fluid, freely falling towards the black hole, the z-component of 
velocity is taken to be u = —y/1 — a 2 . 

Stiff fluid: p = p = p/2. 
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Using these assumptions, we solve the determinant of the coefficients of con- 
stants of Eqs. (14.in - fj4.5j) which results a complex dispersion relation [26J. 
Comparing the real and imaginary parts, two dispersion relations are ob- 
tained. The real part gives a relation of the form 

A 1 (z,u)k + A 2 (z,u) = (4.6) 

which is linear and yields only one value of k. The imaginary part gives a 
dispersion relation quadratic in k, i.e., of the type 

B^^k 2 + B 2 (z,u)k + B 3 (z,u) =0. (4.7) 

Using the values of k, we can calculate the refractive index and its change with 
respect to angular frequency which helps us to study the wave properties. 

Figure 1 shows the graphs of k obtained from Eq. (14. 6j) whereas Figures 2 
and 3 represent values of k obtained from Eq. (j4.7j) . The graph labels A, B 
and C denote the graphs of the wave vector, the refractive index and change 
in the refractive index with respect to angular frequency us respectively. 

Here we are give a brief description of Figure 1 which helps to access the 
results for the other figures. 

In Figure 1, the wave vector admits positive values which indicates that 
the waves move away from the event horizon. The refractive index is greater 
than one throughout the region and increases in a small region 1.2 < z < 
1.8, < uj < 10 with the decrease in z. The change in refractive index with 
respect to angular frequency is positive at random points which shows that 
the waves disperse normally at those points. The dispersion is anomalous at 
rest of the points due to negative values of 

The information obtained from Figures 1-3 is given in the following table. 



Table I: Direction and refractive index of waves 





Direction of waves 


Refractive index (n) 


1 


Move away from the event horizon 


n > 1 and increases in the region 
1.2 < z < 1.8,0 < lu < 10 
with the decrease in z 


2 


Move towards the event horizon 


Same as above 


3 


Same as 1 


Same as above 



In Figures 1,2 and 3, dispersion is normal as well as anomalous at random 
points. 
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5 Rotating Non- Magnetized Hot Plasma 

This section is devoted to rotating non-magnetized hot plasma flow. When 
we consider non-magnetized plasma in rotating background, i.e., B = 0, 
the equations of evolution of magnetic field (13.21) and (13.31) are satisfied. We 
substitute C5 = 0, B = 0, A = and cq = in the Fourier analyzed perturbed 
GRMHD Eqs.f l3T28l) - fl3~3Tl) and obtain 

ci(— p) + c 2 (—p) + c 3 (p + p)[—-i 2 u + (1 + 7 V)tA; 
a a a 

-(1 - 2 7 V)(1 + 7 V)- + 2 7 WV'] + c 4 (p + p) 7 2 [(— 

u ' a 

+iku)V + u{l + 2 1 2 V 2 )V + 2 1 2 u 2 Vu'} = 0, (5.1) 
ci[p 7 2 w{(l + 7 V 2 )V + j 2 Vuu'} + >y 2 Vu(p' + Lkp)} + c 2 [p 7 2 u 

x{(i + rV 2 )\/' + 7 V™'} + rVuQo' + **p)] + c 3 [(p + p) 7 2 

x {(— + fcArafrVii + (1 + 2 7 V)(1 + 2 7 V 2 )V r/ - 7 W 
a 

+2 7 2 (1 + 2 1 2 u 2 )uVu'} + 7 V(1 + 2 7 V)(p' + p')} + c 4 [(p + p) 

{ 7 2 (1 + 7 V 2 )( — + iku) + 7 4 m((1 + 4 7 V 2 W + 4YV'(1 + 

7 V 2 ))} + 7 \(1 + 2 7 V 2 )(p' + p')} = 0, (5.2) 
ci[p 7 2 {a 2 + (1 + 7V) W + 7 VW} + 7 V(p' + 
+c 2 [p 7 2 {a z + (1 + 7 V)mt' + 7 VW} + (p + tfcp) x 

(1 + 7 V)] + c 3 [(p + p) 7 2 {(l + 7V) (— + iku) + u'(l + 7 2 n 2 ) 

a 

x (1 + 4 7 V) + 2u-f 2 (a z + (1 + 2 7 V) VV')} + 2 7 2 m(1 + 7 V) 

x(p + p')] + c 4 [(p + p)7 4 {(— + tku)uV + u 2 V'{\ + 4 7 V 2 ) 

a 

+2V{a z + (1 + 2 7 V)W)} + 2 7 W(p' + p')] = 0, (5.3) 

cip 7 2 [^ + «{a, + ^{VV + uu')}} + c 2 p[— ( 7 2 - 1) 

a a 

+-1 2 u{a z + -i 2 (W + W)}1 + c 3 (p + p) 7 2 [ 7 2 m(^^ + t/cw) 

a 

+(a z + 2 7 W)(1 + 2 7 V) + 7 2 VV\1 + 4 7 V)] + c 4 (p + p) 7 2 

x + tA;it) + 7 2 u\Z'(l + 4 7 V 2 ) + 2 7 2 uV> 2 

+2 7 2 W)] = 0. (5.4) 
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These equations can be used to find dispersion relations. 
5.1 Numerical Solutions 

For rotating plasma, the velocity assumption can be modified as follows. 
• Velocity components: u — V, x and ^-components of velocity are mod- 



ified as u = V = — 




Thus, the Lorentz factor becomes 

7 = 1/Vl -u 2 -V 2 = 1/a. 

This assumption along with the assumptions of time lapse, specific enthalpy, 
density and pressure given in Section 4, satisfy the GRMHD Eqs. (13. 2\i - (13.61) 
for the region 1.4 < z < 10, < uj < 10. The determinant of the coefficients 
of constants in Eq. (l5.1l) - (l5.4p leads to two dispersion relations. The equation 
obtained from the real part of the determinant is quartic, i.e., 

A^zjk 4 + A 2 {z, uj)k 3 + A 3 (z, uj)k 2 + A 4 (z, u))k + A 5 ( z, lu) — (5.5) 

which yields four values of k out of which two are real and other two are 
complex conjugate of each other. The imaginary part gives a cubic equation 
of the form 

B^k 3 + B 2 (z, uj)k 2 + B 3 (z, uj)k + B±(z, uj) = (5.6) 

which leads to three roots in which one is real and other two are complex 
conjugate of each other. 

The real values of k obtained from Eq. (15.51) and (15. 6p along with respective 
refractive index and its change with respect to angular frequency are shown 
in Figures 4, 5 and 6. 
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The results derived from the Figures 4-6 can be represented in the fol- 
lowing tables. 



Table II. Direction and refractive index of waves. 





Direction of waves 


Refractive index (n) 


4 


Move towards the event horizon 


n > 1 and increases in the region 
1A< z < 1.8,0 < uj < 10 
with the decrease in z 


5 


Same as 4 


Same as above 


6 


Same as 4 


n > 1 and increases in the region 
1.4 < z < 1.7,0 < uj < 10 
with the decrease in z 



The regions of normal and anomalous dispersion can be summarized as 

Table III. Regions of dispersion. 





Normal dispersion 


Anomalous dispersion 


4 




1.4 < z < 10,0 < u < 10 


5 


1.4 < z < 10,0 < u < 10 




6 




1.4 < z < 10,0 < w < 10 



6 Rotating Magnetized Hot Plasma 

Here, we assume that the plasma is magnetized and rotating. The velocity 
and magnetic field of fluid are assumed to lie in xz plane. The respective 
Fourier analyzed perturbed GRMHD equations will remain the same as given 
by Eqs. fl3^5D - (l33H in Section 3. 

6.1 Numerical Solutions 

We shall use the same values of time lapse, specific enthalpy, density, pressure, 
x and ^-components of velocity as given in Sections 4 and 5 with the following 
restrictions on the magnetic field. 

• When we take u = V and V F = 0, Eq.([33]) leads to A = 1. 



• B = 
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These restrictions satisfy the GRMHD Eqs.(JS2])-(JSSD for the region 1.4 < 
z < 10, < lu < 10. From Egs- flX^]) and (ET27I) . we have c 5 = 0. Substitut- 
ing C5 = in Eqs. (13.251) and (I3.28p - fl3.3ip . we obtain a 5 x 5 matrix of the 
coefficients of constants. Its determinant leads to two dispersion relations. 
The real parts give 

Ai(z)£; 4 + A 2 (z, uj)k 3 + A 3 (z, u)k 2 + A 4 (z, u)k + A 5 (z, u) = (6.1) 

yielding four values of k out of which two are real and interesting shown in 
Figures 7 and 8. The imaginary part leads to 

Bi(z)k 5 + B 2 (z, to)k A + B 3 (z, cu)k 3 + B A (z, cu)k 2 + B 5 (z, u)k 
+B 6 (z,u) = 0. (6.2) 

This represents fifth order equation giving five roots in which three are real. 
These solutions are represented by Figures 9, 10 and 11. 

Results related to the direction of waves and the refractive index are given 
below. 



Table IV. Direction and refractive index of waves. 





Direction of waves 


Refractive index (n) 


7 


Moving away from the event horizon 


n > 1 and increases in the region 
IA < z < 1.8,0 < u < 10 
with the decrease in z 


8 


Same as 7 


Same as above 


9 


Moving towards the event horizon 


n > 1 and increases in the region 
1.4 < z < 1.7,0 < lu < 10 
with the decrease in z 


10 


Same as 9 


Same as above 


11 


Same as 9 


Same as above 
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The information about the regions of normal and anomalous dispersion 
obtained from Figures 7-11 is given in the following table. 



Table V. Regions of dispersion. 



Figure 



Normal dispersion 



Anomalous dispersion 



3 < z < 10, 1 < uj < 10 
2 < z < 10,2 < uj < 2.4 
2 < z < 10,3.41 < uj < 3.6 
2 < z < 10,4.2 < uj < 4.6 
2.3 < z < 10,5.4 < uj < 5.62 
2<z< 10,6.38 < < 6.6 
2<2<10,7<cc>< 7.35 
2 < ^ < 10,8 < uj < 8.2 
2 < ^ < 10,9 < uj < 9.2 
2 < z < 10,9.9 < uj < 10 



1.4 < z < 10,0 < w < 10 



1.8 < z < 10,2.8 < uj < 3 
1.9 < z < 10,3.6 < < 3.8 
2.3 < z < 10,4.68 < uj < 5 
1.6 < z < 10,5.65 < < 6 

1.8 < z < 10,6.6 < uj < 6.8 

1.9 < z < 10,7.38 < < 7.6 
2 < z < 10,8.3 < uj < 8.97 

1.8 < z < 10,9.58 < < 9.84 



2.1 < z < 10, 1.7 < uj < 2.3 
2<z< 10,4.21 < < 4.6 
2.3 < z < 10,5.1 < < 5.4 
2 < z < 10,5.67 < u < 6.1 
2 < z < 10,7.21 < < 7.6 
2 < z < 10,8.01 < u < 8.2 
2 < z < 10,8.4 < < 8.6 
2<z< 10,9.41 < uj < 9.65 



2 < z < 10,2.5 < < 2.8 
2<z< 10,4.61 < uj < 5 
2 < z < 10,5.41 < < 5.65 
1.8 < z < 10,6.14 < uj < 6.4 
1.8 < z < 10,6.9 < uj < 7.2 

1.7 < z < 10,8.6 < < 8.88 

1.8 < z < 10,9.69 < uj < 9.92 



10 



1.7 < z < 10,1.21 < uj < 1.4 
2 < ^ < 10,2.6 < < 2.77 
1.9 < z < 10,3.4 < < 3.6 

2 < ^ < 10,4 < uj < 4.2 
2 < z < 10,4.61 < < 4.97 
2 < ^ < 10,5.4 < uj < 5.69 
2 < z < 10,6.9 < uj < 7.2 
2 < z < 10,7.6 < uj < 7.8 
2<z< 10,8.01 < uj < 8.35 
2 < z < 10,9.2 < < 9.4 



1.7 <z < 10,0.4 < uj < 0.57 
2 < z < 10,3.65 < < 3.85 

2 < z < 10,5.7 < uj < 6.1 
2 < z < 10,6.44 < uj < 6.8 
2 < z < 10,7.3 < < 7.6 
2 < -2 < 10, 7.8 < uj < 8 

1.8 < z < 10,8.37 < uj < 8.6 
2.2 < z < 10,9.5 < < 9.81 



11 
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7 Summary 



In this paper, we find the wave properties of hot plasma in Schwarzschild 
magnetosphere. For this purpose, we have derived the 3+1 GRMHD equa- 
tions for the scenario. Their component and Fourier analyzed forms are 
formulated using the specific assumptions. Dispersion relations are calcu- 
lated for the non-rotating, rotating non-magnetized and rotating magnetized 
plasmas. 

For hot plasma living in non-rotating plasma, our assumptions satisfy the 
3+1 GRMHD equations in the region 1.2 < z < 10. In Figure 1, we have 
found that the waves are directed away from the event horizon while Figure 2 
shows that the waves move towards the event horizon. In Figure 3, the waves 
are directed away from the event horizon. In all these figures, the dispersion 
is found to be normal and anomalous randomly. 

For the rotating non-magnetized plasma, the assumed parameters satisfy 
the 3+1 GRMHD equations in the region 1.4 < z < 10. All the figures 
indicate that the waves are directed towards the event horizon. In Figure 
4, dispersion is anomalous in the whole region while Figure 5 indicates that 
dispersion is normal throughout the region. In Figure 6, dispersion is anoma- 
lous. 

For the rotating magnetized plasma, our assumptions satisfy the 3+1 
GRMHD equations in the region 1.4 < z < 10. In Figure 7, the dispersion 
is anomalous and waves are directed away from the event horizon. Figure 
8 shows that the waves move away from the event horizon and dispersion 
is normal in most of the region. Figures 9, 10 and 11 admit random points 
of normal dispersion. In these figures, the waves move towards the event 
horizon. 

The summary of results and comparison with previous literature is given 
as follows. 

We compare our results with the previous work on isothermal plasma [23J . 
We have used here the variable specific enthalpy while in previous work this is 
constant. Here propagation vector admits both positive and negative values 
which shows that the waves can move away and towards the event horizon. 
In the previous work, this always takes negative values which indicates that 
the waves move towards the event horizon. This comparison shows that 
variation in specific enthalpy effects the direction of waves. The refractive 
index is always greater than one in both cases. Dispersion is normal and 
anomalous at random points as shown by figures in both works. Here the 
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waves move away from the event horizon and disperse normally in most of 
the region of figure 8. This shows that there is a chance to obtain information 
and energy from the magnetosphere. In previous work, waves always move 
towards the event horizon which means that no information can be extracted 
from black hole whether dispersion is normal or anomalous. The comparison 
of the results of cold and hot plasma shows that gas pressure grows the 
normal dispersion. 



Table VI. Comparison of results. 



Results 


Previous work 


Recent work 


Direction of waves 


Towards the event horizon 


Away from the event 
horizon 


Dispersion 


Normal at random points 


Normal in most of the 
region in Figure 8 


Conclusion 


No information can be 
extracted 


A chance to obtain 
information from 
magnetosphere 



It is concluded that the hot plasma waves in the Schwarzschild magneto- 
sphere can have an escape towards the outer end of the magnetosphere. 

Acknowledgment: We would like to thank Dr. Umber Skeikh for the 
fruitful discussions during this work. 



Appendix A 

This appendix contains the GRMHD equations for the general line element 
and the Schwarzschild planar analogue. The set of Maxwell equations and 
conservation laws under the influence of gravitational field are collectively 
known as GRMHD equations. The hot plasma model includes the mass, 
momentum and energy conservation laws. The GRMHD equations for the 
general line element, given by Eq. (12. 1 j) . are [23] 

^ + -(B.V)/3 + #B = -Vx (qVxB), (Al) 
ar a a 

V.B = 0, (A2) 
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DP ° + Po7 2 V.^ + * I Z± + V.(aV - (3) } = 0, (A3) 



B 2 \ i 1 1 /)yi 

Po/V + i- 7« + Po^ViYj - —BiBj 



+Pol 2 V^ ~ (^7« - Vf fe V fe = - Po 7 2 pR 

- (^7^'} + ^(V x B),V.(V x B) 

(aV x B) - (B.V)/3] + (V x B).V/3})i, (A4) 
^(tVpo) - | + ^((7 2 PoP - P ) + ^((Vx B) 2 + B 2 )) 
+ ^{7W^^ +p 7 y ' - ^-((V x B)*(V x B)^' 

Za 47T 

+ ^((V x B) 2 + B 2 )yi}£ tlij + PPo7 2 V.a 
+/i P0 7 4 V.(V.V)V - -U.((V x B) x B) - -L(V x 

Z7T 47T 

(V x B)).B - - ( m7 2 V.(V.V)/5 + p(V./?) - -i-B.(B.V)/? 
a [ 47r 

(V x B).{(V x B).V}/3 +t^((V x B) 2 + B 2 )(V./?)1 

57T J 



1 

47T 



1 

47T 



(V x B).(a x B) + 2(V x B).-^-(V x B) 

dr 



— (V x B).{(V x B).V}/3 + 0(V x B).(V x B) 

a 

dr 



= 0. (A5) 



For the Schwarzschild planar analogue, f3, 9 and £t r Yij vanish, the perfect 
GRMHD equations take the following form 

_ = Vx(aVxB), (A6) 
V.B = 0, (A7) 
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do dV 

-g + (aV.V)po + Po7 2 V.^ + p 7 2 V.(aV.V)V 

+p V.(«V) = 0, (A8) 

{(/W + 5-)^. + p p 7 V^ - i-s^-X— + v.v)v* 

- ( 4^' - l B ^ yj > k yk + ^ 2l ^f + (V - VM 

= -po/iT 2 ^ - P„ +^(V x B),V.(V x B) - -i^(«B) 2 n 

47T 87rcr 
+-L(aB i ), j -l[Bx{Vx(Vx(aVx B))}],, (A9) 

7 2 -|(m) - + PoP7 2 (V.a) + p p 7 4 V.(V.V)V 

a at a at 

_i_ a .{(V x B) x B} - -L(V x (V x B)).B 
+ — <^ 2-(V x B). — (V x B) + (V x B).(a x B) 

47T ^ OL at 

+ B.if}=0. (MO) 
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Surface Plot of of Wave Vector 




S ux i ice Plot of Refract ive Index 



Surface Plot of D[n 




Figure 1: The waves are directed away from the event horizon. The dispersion 
is found to be normal and anomalous randomly 



ABC 




Figure 2: The waves moves towards the event horizon. The dispersion is 
normal at random points 



ABC 




Figure 3: The waves move away from the event horizon. The dispersion is 
found to be normal and anomalous at random points 
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ABC 




Figure 4: The waves are directed towards the the event horizon. The region 
shows anomalous dispersion 




Figure 5: The waves are directed towards the the event horizon. The disper- 
sion is normal 



Surface Plot of of Wave Vector 




Figure 6: The waves are directed towards the the event horizon. The whole 
region admits anomalous dispersion 
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Figure 7: The waves move away from the event horizon. The whole region 
admits the anomalous dispersion 
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Figure 8: The waves are directed away from the event horizon. The dispersion 
is normal in most of the region 
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Surface Plot of of Wave Vector 
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Figure 9: The waves are directed towards the event horizon. Normal disper- 
sion is found at random points 
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ABC 




Figure 10: The waves move towards the event horizon. The whole region 
admits the random points of normal and anomalous dispersion 




Figure 11: The waves are directed towards the the event horizon. The dis- 
persion is normal at random points 
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